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ABSTRACT 


This thesis studies the problem of finding the 
actual size of a Welch-Aspin test as applied to the 
solution of the well-known Behrens-Fisher problem. 
Methods of finding the actual size, due to Wang, and 
Mehta and Srinivasan, are discussed. Modifications to 
Wang's method are attempted, and the resulting methods 
are used to find the actual size for a larger number of 
parameter values than used in previous studies. 
Numerical results are given for each method. These 
results are compared, not only with each other, but 
also (where applicable) with those results obtained in 
earlier studies. It is seen that the modified method of 
Wang tends to produce a size of test in closer accordance 


with an assigned size of test a. 
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CHAPTER 1 
INTRODUCTION 
Lat The Problem 


im thestield of Statistical data analysis, 1 as 
often necessary to compare two populations with respect to 
pecentain, Characteristic. For exampte, agiiculturalists 
May want to compare the yields of wheat from two different 
farms. Examples such as this form the basis of the Two 
Means Problem, which is defined as that of comparing 


the means of two normal populations. 


1 is the mean, and Say the sample 


variance, of a sample of size nj, taken from a normal popu- 


Suppose that x 


lation with true mean Wy and true variance oie and that 


Kor Bas Ny, Uns and of are the corresponding values from 


a second normal population. 


If the exact values of the population variances 


ae and os are known, then the problem of testing the 
hypothesis Hy = U5 can be solved using the Normal dis- 
tribution. If the exact values are not known, but the 


GabLOws a= o5*/ oe is known, then this problem is solved 
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which has the t-distribution with n, +n 


1 ~ 2 degrees of 


2 
freedom if Wy = Uo: 


If, however, the value of 6 is unknown, then the 
above d statistic cannot be calculated: the solving of the 
two-means problem when @ is indeed unknown, is known as the 


Behrens-Fisher problem. 
le2” ReEView Of the Literature 


The solving of the Behrens-Fisher problem has been 
the subject of a great deal of research. It was first 
considered by Behrens (1929), who suggested that the dis- 
tribution of the difference between two means could be 
expressed in terms of observations in the samples from the 
two normal populations. Fisher (1936) extended this as 
the correct "exact" solution and Sukhatme (1938) drew up 
tables for the distribution of the Behrens-Fisher statistic. 
Bartlett (1936, 1939), among others, criticized this 
solution by saying that the probability of rejecting the 


hypothesis of the equality of the two means (using Fisher's 
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5 
method) was generally less than the assigned size of the 
test, and tried to solve the problem using the theory of 
confidence intervals. Welch (1938) briefly mentioned this 
method, and Scheffé (1943) wrote up the full solution: 
this solution was based on the t-distribution, and was 


proven CO be the best. solution of its “fype available. 


‘Welch (1947) developed a solution based on a 


series approximation for a critical value Y ces such 


that 
> = 
BY ane Vv, te) } OL 
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Aspin (1948), as well as Trickett, Welch and James 
(1956), tabulated critical values Vibe) fOomPtnius solucion for 
certain sets of values a, Ny n5 ana ce sOther solutions 
have been proposed by Wald (1955), Banerjee (1961) and 
Pagurova (1968). Scheffe (1970) provided a recent paper 


that reviewed many of the above solutions. 
devo Objectives of the Research 


In the paper of Scheffé (1970), the question is 
raised as to the actual size of the test for Welch's solu- 
tion. Welch, himself, in an appendix to Aspin (1949) 
found the actual size for a few values of the parameters 
= ia) shy) eh ot OR Abe 

9 A 
specified method. Mehta and Sr 


tos 
inivasan (1970), while 


, 3 2, Ly using san an— 
comparing the actual size of the test for various solutions, 
used a method of triple integration developed by Golhar 
(1964). Wang (1971) found the actual size of the Welch 

test for selected sets of parameter values by using a method 
involving numerical quadrature and the t-distribution. How- 
ever, all of these writers presented results for only a few 
sets of a, ny and No- Therefore, it is the objective of 
this research to investigate these methods for the purpose 


o£ determining the actual size of the test for a greater 


number of sets of parameter values. 
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CHAPTER. 1I1 
WANG'S METHOD 


In the paper of Wang (1971), it vis stated that for 


apgiven test withestatistic w =) (x = 1) > cand critical value 
S 


Be (c) at an assigned size of test o, the actual size of 


the test may be written as 


P (C78) = cate =e Cy (c) 0 ) 
S 
where 
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Multiplying both sides of the inequality by the 
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C, fe) = Welch-Aspin critical value 
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Using the hypothesis that Hy = Us the left- 
hand side of the inequality in the above expression is a 
Student's-t variable with fy - f. degrees of freedom, while 
the right hand side is a function of the random variable 


r and is distributed independently of t. 


Wang defines a Beta (Cte a variable z as 
2 2 2 
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This is equivalent to saying 


2 
z f, (£541) oF 


‘spe io (Jae) 
(1-2) f, (£,+1) o5 


TiewlLdieenandsside: of the yinequal Leyei(2.4.) <s thus 
labelled as a function of the Beta variable z, say h(z), and 
the complete expression (2.1) is written as an expectation 
over the Beta variable z 


raf ey 


ey E, | Pr [te 4s, 2 h(z) | 0,2] } 
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where one {h(z)] is the upper tail of the t-integral (to 


thescigne cr. therpoink (Zz) ) with f,+f, degrees of freedom. 


The value of the actual size of the test can then 
be calculated by quadrature using Simpson's Rule in the 


LOLLOwing, forme 


2n 
eae jem es" 
P (C,,8) = E(O¢ +e, b nla) = 30 Cee eo 
2n = number of intervals to be used in Simpson's Rule. 
z. = — (note that 02,8 1 as Ze is a Beta variable) 
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An algorithm to evaluate this expression is given 


as follows: 


1. Choose the values of a, Ny, Ny and 8 to be 


tested (there will be nine values of SoS tested for 


each set of Ny, Ny and a, namely l, 
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2. Fit a fourth-degree polynomial to the avail- 
able 11 values in the particular row of a Welch-Aspin 
table of critical values corresponding to the choice of 


a, Ny and Nn. (this polynomial will be used to find the 


value of Cle), the Welch critical value, necessary to 


calculate a particular value of h(z;)). 
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3. Decide how many intervals will be used, and 


Cari this value ei. 


4. For each value of Zs (A9=20 EO: 2h) carculare 
the value of h(z,) by transforming Zs to rs using (2.2) 
and making the appropriate substitutions in the right- 
hand side of the inequality in (2.1). Use this value of 
h(z, ) to determine the value of the t-integral to the 


TIght of the point h(z; ) namely, QF +e, Ch (2, )] : 


5b. gOnce all of the (2nt)) values or Or, +£,! h(z,)] 


have been determined, use Simpson's Rule in the form given 


earlier to evaluate E Q flezae) | , which as the 
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Gesired estimate of the actual size of the test. 
Two things should be mentioned here: 


a), The fitting of the polynomial to the 11 points 
in a Welch-Aspin table of critical values corresponding to 


particular values of Ny, Ny and a can be done using any 


accurate curve-fitting method. The polynomial should be 
fittedavusing the values of c indicated in the tables (c= 
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CHAPTER Iii 


MEHTA AND 


SRINIVASAN'S METHOD 


in the paper,of Mehta and Srinivasan (1970), the 
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an and variance of a second normal 
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general form. 
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and X,1 X5, are the means and Sy ; S.5 the variances of 
samples of size ny and no taken. from the two normal popu- 
lations. 
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expressed as (for an assigned size of test a) 
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The power of the test is defined as the probability 
of rejecting Ho when 2t.i1s false; or in other words, Prob 
PSS Vi fe) f or 1- g(n,9). This power is equal to the 
actual size of the test when n =.0. Therefore, the actual 
Sizes0t ehe test 15 diven by the value of l= 60,0). fhe 
calculation of this value necessitates the evaluation of 


the triple integral (3.2) with p replaced by 0. 


The inner integral in (3.2) is simply the standard 
normal cumu lativesarstribution function © (W(x,;y) ), and can 


be evaluated by using the relationship: 
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c + erf (W(x,y)) ) dxdy . Ca 
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This resulting double integral can be evaluated using 
Simpson's Rule, by making the following transformations on 
x and y in order to make the limits finite, and using the 


FORTRAN SLunCeLon ERE scOlvevValuate, the Crron f£unccion: 
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From examining this integral, we can see that it 
issundetaned at s—l, or s =I) vand close examination of the 


formula for W ( Ss Sa will Show that it is undefined 
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at s=0. 


Therefore, a small value of € (Say, 1077) is chosen, 


and the value of the integral is found using the limits 
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The Simpson's rule form for the above integral is 
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h = mi (the length of each interval on the 
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k = ae (the length of each interval on the 
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M = number of points along each axis 
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= weighting coefficients which can be found 


TABLE OF WEIGHTED COEFFICIENTS FOR TWO-DIMENSIONAL 
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CHAPTER 1V 
PROGRAM AND RESULTS 
4.1 Wang's Method 


For Wang's method, a program was written in 
FORTRAN, uSing the procedure as outlined in the 
algoricnin given an Chapter ff (a flowchart ton this 
program is given in Fig. 4.1). As Wang's method involves 
direct use of the Welch-Aspin tabulated critical values, 
it was decided to find the actual size of the test as 
such ior the same values of a (and all finite pairs of 


sample sizes Nyy n,) used in setting up the Welch-Aspin 
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1+6 
The results are given in Table 4.1. Investigation of 


these results shows that in all cases, the tabulated 
size (to three decimal places) is equal to the assigned 


Sige Of FEST "a. 
4,2 Mehta and Srinivasan's Method 


For the Mehta and Srinivasan method, a program 
was written in FORTRAN using the procedure as outlined in 


Chapter-111l”* (a flowchart-for this.program isigiven in 
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START 


Read in a row of a Welch-Aspin Table, 
and values of fj, f2. 


Fit a 4th degree polynomial to Welch-Aspin 
see? Gleheel Wisishiove; WreWlWelsy tone C(O, oP, 056-51) - 


Pick a value for M(the number of 
points to be used in Simpson's Rule). 


Select ae Oo 


Use the formula of Wang's 
to calculate the value h(z,). 


Use the Series Approximation of Hill 
to calculate Q, LHe ye 
itf2 


NO 


Use Simpson's Rule to find 
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FIG. 4.1 FLOWCHART FOR WANG'S METHOD 
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TABLE OF CALCULATED ACTUAL SIZES 
USING WANG'S PROCEDURE 
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TABLE 4.1 (continued) 


TABLE OF CALCULATED ACTUAL SIZES 
USING WANG'S PROCEDURE 
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TABLE 4.1 (continued) 


TABU OF CALCULATED ACTUAL SIZES 
USING WANG'S PROCEDURE 
(a = .01) 
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TABLE 4.1 (continued) 
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USING WANG'S PROCEDURE 
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Bagune, 422) 4). 3). galt was; decided to ifindrthe actual size of 
thestest. using tthewsame valveseof wi(andeall finite pairs 
of sample sizes Nyy n.) used in setting up the Welch-Aspin 
San les, Biome Mes. Cie (of lc000 1 F 00d €.01l a ws, 1, 110), 
100, 1000, 10000) which Mehta and Srinivasan had used pre- 
viously. In addition, values were computed for the 


pawureu Ear weasem0n— 05, Wl. = i. = 7 ,and 0 values of (5, 


1 2 
a = a 2" o i 3! +) in order to be able to compare 
the results of Mehta and Srinivasan's method, Wang's 
method, and Welch's unspecified method. The results for 
the general testing of Mehta and Srinivasan's method are 
given in Table 4.3, and the results for the special case 
are given in Table 4.4. Investigation of these results 
shows that in all but a few cases, the tabulated size (to 
three decimal places) is equal to the assigned size of 


test a. In all cases, the actual size to three decimal 


places is within +0.001 of the value of a. 


It should be noted that, due to the limitations of 
the IBM 360/67 computer being used, results were not ob- 
tained for any set of parameters where fy was equal to 30. 
This is because the upper limit of s used. (.9999999) was 
equivalent to x = 9999999, and the calculation of the term 

Clee t) . 
pay as for these particular values of x and ft) caused 


exponent overflow. Lowering of the upper limit of s to 


,99999 was tried) but this produced a highly inaccurate 
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Tesvitetor thepactial sizemot the test. For simi Lar 
reasons, results for fr, = 30 were also unobtainable. 
Therefore, in order to get an idea of how accurate the 
method actually is for values of ft, (or f.) = 30, 
another, more time consuming, version of the Mehta and 


Srinivasan method was implemented, as described in the 


following, algorithm. 


1. Decide on values of fo. fo, and a to be used. 
2 


g 
2. Choose a value of § = a , and M (number 
2 
Of points) . 


3. Set NN equal to 1 and AVAL = 100.0. 


4. Set the upper limit of the integral equal to 


2 e 
5. Use Simpson's Rule in two dimensions, to 
evaluate. 
NN NN plvcga’ tb labo 
2 2 oS leone a 
SUM = ie, e x y (l+erft (W(x,y)) dxdy. 
€ 0 


using M points in each dimension. 


ae If the absolute value of SUM-AVAL is less than 
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ey then go to step 11.) Else 
7. Replace AVAL by the value of SUM. 
8. Increase NN by l. 
9. Change the value of M to 2(M-1) + 1. 
10. ~GO°to step 42 
ll. Set the actual size of the test equal to 
SUM 
ae z 


iL 
20s yl ts) 


12. If all values of theta have not been tested, 


then go to step 2.~Else 


eo LOD. 


As computer implementation of this version of the 
Mehta and Srinivasan procedure is extremely time consuming 
and costly, it was decided to use the above algorithm for 
only one set of data (fy = 10, f. =—'§30,.°0 =...01)..) Results 


were as follows: 
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TABLE 4.2 
CALCULATED SIZE OF TEST 


(MEHTA AND SRINIVASAN PROCEDURE) 
ST ns eer — re 0 es eee — oie 5 6) 


actual 
Stems OrOO mw. 02 00 OL 0. 0200 O02 Ol O21 O01 Od Ones LO od: 


It can be seen that these results are also in agreement 


(to three decimal places) with the value of the assigned 
Size of test a. The author believes that similarly 
accurate results will occur for other sets of parameter 


values with one or both of n equalling 30. 


4.3 Comparison of Results 


Table 4.4 shows a comparison between the actual 


sizes of the test for eh at 7) Oo) e005, sand. Va kuesm@or 
1a Cd re car Ne, 
6 of (= — =’ =), as calculated by Welch (1949), 


2, 
ORO cea), Komal 
Wang (1970), and the present writer (using both Wang's 


and Mehta and Srinivasan's basic methods). 


The results of Welch and the present writer 


START 


Determine starting value of M (no. 
of points in each direction to 


be used for Simpson's Rule). 


Determine value of assigned size 
Ope INS Cl, chore) seslinvel seas Che (Sc. 


Read values of Nl and N2. 
Calculate le (Nii) ae le(NiZ a) 
2 2 


PrckevValuewormc. 
Calculate C= Ay 


A= 1/N; 


Set AVAL 


= 100, 
MM = S 


tarting value of M. 


Use Simpson's Rule to evaluate transformed 
Mes Ie 
integral i f *** using MM points in each direction. 


€ 
Calbliaares SUM: 


— NO 
AVAL=SUM | S(AVAL-SUM) & & 


MM = ((MM-1)*2)+1 ies Seen ce 


| YES 


et actual size of test =1l - SUM 
27 (NL-1) Ens) 
2 2 


ANnGs pis liEmO Wis 


<——All values of @ tested? 


— pid Sets Of Ni, NZ teszed?s = 
Se 


— All values of a tested? NO 
FINISH 


FIG. 4.2 FLOWCHART FOR MEHTA AND SRINIVASAN'S METHOD 
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Calculate Simpson's 
Rule coefficient W(I,J). 


Calculate c = al 
ay ¢ eye GL e Low 
Te) {f)/£2)() 


Calculate Vv, (©) using Welch Series Expansion. 


Calculate W(x,y) = Vv, (e) eure 2y-Gier)., 
fy nD 
Calculate ERRORF = ERF ,W(x,y). 
(et) 
3a) ie 


2 
aR Vers a ee poe 
VAL =e x? y * (1+ERRORF'). 


(1-s)2(1-s')2 


M-1 M-1 
Leste) aslo gnc 


NO 


SEas + ae 


NO 


EXIT 


FIG. 4.3 FLOWCHART FOR SIMPSON'S RULE 
(AS USED WITH MEHTA AND SRINIVASAN'S METHOD) 
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TABLE 4.3 


TABLE OF CALCULATED ACTUAL SIZES 
USING MEHTA AND SRINIVASAN'S PROCEDURE 
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TABLE 4.3 (continued) 


TABLE OF CALCULATED ACTUAL SIZES 
USING MEHTA AND SRINIVASAN'S PROCEDURE 


(a = 12025) 


20252 
2025. 
Oo 
70250 
ALC PASH 


e0252 
AA PsN 
AO SyAP 
OPS 
ORAS 


Oia 
402 Srl 
AUPAS 
70250 
702 52 


«252 
OZ od 
eae el 
ue Od 
eo 


ee Onc 
nO Se 
MOP AS, 
L252 
Uo 


20252 


FU Zou 
OW Ase 
OPASY 
50252 


teed 
ica 
0250 
O20 
~0252 


ia 
eo ae 
a0 PLONE 
~0250 
e022 


W252 
Oras 
OMe al 
2 U2 oi 
Oe 


0203 
Oey 
U2 2 
Oe oe 
LUZ 53 


act Los rae _ 7 = _ ; 7 
| cca ts aie 
. 7 : a 7 i a 7 . 
- 7 : i: : 
> : (bauats: oe) Yh 4A? : 


7 QRkiz MAUTSS GUPRIATD TO Bae | 
qHyaro &° Terre oHA «Tee oat 


33 


TABLE 4.3 (continued) 


TABLE OF CALCULATED ACTUAL SIZES 
USING MEHTA AND SRINIVASAN'S PROCEDURE 
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TABLE 4.3 (continued) 


TABLE OF CALCULATED ACTUAL SIZES 
USING MEHTA AND SRINIVASAN'S PROCEDURE 
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TABLE 4.4 


ACTUAL SIZES OF TESTS: COMPARISON OF RESULTS 


Welch Present 


(1949) (1970) Applied Applied 
Wang's Mehta's 


(Wang's method) are in close agreement (both with each other 
and with the assigned value of a). However, there is a 
Slight discrepancy between these results and those of Wang 
and the present writer (Mehta and Srinivasan's method). 

The possible inaccuracy of Wang's results (1971) could be 
due to the fact that she may not have used as accurate an 
approximation (such as that of Hill (1970)) for the evalua- 


tion of the upper tail of the t integral. Any inaccuracies 
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of the results of the present author (Mehta and Srinivasan's 
method) are undoubtedly due to the way in which this proce- 


dure was implemented on the computer. 
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CHAPTER V 
SUMMARY AND CONCLUSIONS 


The methods of Wang and Mehta and Srinivasan for 
finding the actual size of a Welch-Aspin test as applied 
to the Behrens-Fisher problem, have been investigated in 
detail. Computer calculations of the actual size have 
been done for a large number of sets of parameters a, Nyy 
No and 8. Of the two methods used, that of Wang was 
easier and more economical to implement on the computer 
for the following reason. Each step in Wang's method is 
clearly defined, and, in many cases, extremely accurate 
computer approximations exist to simplify the calculations 
involved. Good accurate computer approximations for the 
majority of steps in the Mehta and Srinivasan procedure do 
nGteexien —(thiusitis particularly true for sttheccalculating, of 
the double integral specified in the method): this means that 
long laborious calculations must be used. It is also true 
thatthe quantucy ie( 0, G)rcalculated an the Mehta 
Srinivasan method covers a larger area than the quantity 
Es hes [h(z) J) calculated in Wang's method. This is 
because the value of 8€(0,8) must be subtracted from l 


torobtainethe sactual size of the test, which«E—(Q heneAay a 
Zi fy + f. 


) 


tot mses intae bre code BRS eMeW Ro atevizon sf? 
boliqus 26, 4253 ntudA-efen 6 te enis Ieqaas em 
nl. hesapigeever nome overt te! doa 26 42i4-enendell a 

Syvex stre’ leydoe octd td’ webiinivateo Tatuqmne> 


_ 


.pf xe BIetgce hg ic aser Io tedauen sete! 5 38 ome ns 

zie DAEW I> 46490,b980 Shoutsem gw, end 30 bag 

TeziGmen sii), fo dtansiqut ot. Tsaimonone rem bee. 2 

el bodjean &° pre ot Gata sont .dosne2 siiwokaek ait” 
sjatioqs. vVieretsko »2ocen vném at. . BRe boeeeee 
arma ioteo “alt PeLqmin o¢ falas avo bssiizovoRh 

siis: 10! » eooifamlineaga® séfierss adeasess Sacer 

ob awibasere neantt arise fice gtioM octal eystte Jo 

Sec ortgsidstss) sitd. 30) sss vinelandii el erdaj ery , 
Jaceta \engsm eig@t >, [Oodtom sft nt bstiistqe, fetpasnk 7 
suta cela (eit 41. \Sse0 aa Jeon anetretlostaoo tebe 
“Saha vaca nz Setelnotiso (8:1 3)-5 ytiaaaap wit 
espana 38 bila deze septa! = carver hodses’ 


38 
represents directly. There is a chance, therefore, that 
the calculating of 6(0,@) without accurate computer ap- 
proximations may produce a less accurate answer than the 
Caleowlating? of ache al Perhaps because of this 
problem, the method of Wang generally produces answers 
in closer agreement with an assigned size of test a. It 
should be noted, however, that Mehta and Srinivasan's 
method does have the one advantage that it can be applied 
for any finite set of parameter values, while the use of 
Wang's method, as Bote de is restricted to those finite 
and a used in the calculating of the 


values of n n 


ie ieee 
Welch-Aspin tables (Mehta and Srinivasan themselves quoted 


results for fy = f., = 3, ft, = f., = 19, and ti See ee 


19 - values not used by Aspin, Trickett, Welch or James 


in setting up their tables). 


A recent paper of Golhar (1972) has just come to 
the present author's attention. In this paper, he suggests 
evaluating the Mehta and Srinivasan triple integral by 
uSing the error function approximation for the inner in- 
tegral, and Gauss-Laguerre quadrature for the resulting 
double integral. Golhar gives results for the actual size 
of the test using this procedure (for a few select sets of 
if and fo), that appear to be in very close accordance with 


i 


an aSSigned size of test a. The use of Gauss-Laguerre 


aioe 
ae 


ora ANAS swans ss@zuane BeSh, a lsbpmORG en 
sius to yeeunved aquired =. tl teh Ah yg to 

auswaris segubeta yhiuvsnap enel 2 bonsem ote 
2% .o Se64 Io ose \bemoiees, He Gr aw menage: 

=" Nessv inte?! bas catat som \navatod ,badea ad 
betiqus mi Ws> JL Jed! eAesnsvGa BAO SP oved aed 
a> sey od¢ Slidw ecules ‘soeparey YO tee etn 


ejini¢ eo8 of Sagoisdger 2f \bonlisue se bere et 


: a 
4/9 it cnissaluoltss se? ni heen a BAS, A +h te 


boduep 2avisetieds neesviaize put 4*de) =ssidss sige 
' 
ad ve = a Bie ef - yt _ rt i. = gt _ Fas 7ay 


Poteu th) doled \Sialss wn. . qeqeA wa aad tun aeulew = @F 
a 
, (eeines ee Gu wines 


ot? amon daub dad (STEL) 260050 3¢ tegeg dasmed A ‘| 
SIZapplt of 4 fSqs9. 210) ni opokineatie sdetios 36 
qe Tsapattt ele rs). ndeewindee Bas eovedl ants 
"ah thet ahi cot hollsiisosye, modboew? ToT ey: 
gasifiers aes ‘xOd emntetbenp efi svysd+seNeo baw. 
ests lmvd9R uit TOs arbmst Gevte aHltes -Leypesne 
Mr ease apblge wis ® gol) ‘eyvhebsay Sint piew. aan 
seaeicas thin ps Aine st 34 sat ad oe ap ae 


59 


quadrature appears to eliminate many of the problems en- 
countered by the present author in implementing the Mehta 
and Srinivasan method on the computer. Therefore, it is 
recommended that this modification to the Mehta and Sri- 
nivasan method be used in any future investigations of the 
Welch-Aspin test. An attempt should be made at the same 
time to generalize Wang's method for all possible sets of 


finite parameter values. 
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APPENDIX A 
SERIES APPROXIMATION FOR STUDENTS! T 


A procedure for finding the value of the upper 
tari, Or Sludents t+ distribution hasbeen given by 11 
(1970). This procedure involves four different approxi- 
mations - the approximation to be used depends on the 


values of t and n. They are given as follows: 


I For Varge n(n s20)4) or noninteger 1, for all 
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is given by the value of 


1 - ers (=) 
ee Nee 
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2. For small n>l1 and moderate t, with n being 


odd, the upper tail of the t-integral is given by the 


value of 
1 y ge {-- a | n-3 \... 
= | arctan iy) + 5 1 + 3b eae 1 + eS 
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3. For small n>1 and moderate t, with n being 
even, the Upper tail of the t-integral 1s given by the 


value of 
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4. For large t, the series 


s = €(n) we 4 4 eee + ew. ehete 
with 

Gn) Fe= (ntl) /2)/ tar Pin72))" 
and 

w = Bee Sus a , 


is summed until a negligible term occurs. The upper tail 


of the t-integral is then given by the value a8. 


The flowchart on the following page describes 
the computer program that is used to implement this pro- 


cedure. 
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Set s = X+(X243Xx) /B /nT (N/2) 
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APPENDIX B 


A DESCRIPTION OF THE WELCH-ASPIN CRITICAL VALUE 
SERIES EXPANSION 


Welch (1947) derived the series expansion for 
Nee) by first considering the general problem concerning 
N, a population parameter, estimated by an observed 


quantity y which is normally distributed with variance 


where hs are known positive numbers and ae are unknown 
variances. Suppose that ow are the estimated variances 
based on f; degrees of freedom, so that the sampling 

GUS tribGEevon gor Ba is chi-square (distributed independent- 


ly of each other and of y) which is 
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The probability statements about y were made similar in 
character to those which W. S. Gosset derived for the mean 


Of sa single samplie of neobservations ('Student," 1908). 
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Welch sought a quantity h to satisfy the equation 
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it being understood that the exponential is to be expanded 
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Welch then proceeds to develop a series solution by 
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where the exponential is to be expanded in powers of D, 


and that these powers are to be interpreted so that 
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The initial approximation is the large sample 


normal approximation 
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Aspin (1948) extended this series to terms of order 
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> 


A es + = 

2 7 
free 
iM oe 


ww 
a 
ry 


b 4 
_ rz 
— i 
; (34h) 
- r* 7 ) 
r 4 
e A c 
[ ae | sy) 
i oP = eh ie 
yi 


Ssirtc,, 16 (anges “ae ais benmedits (1 
bier .NIZBT) 5 oh late s45h 
ee hi 
7 7 a - 
oe wo 


56 
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(t.e. upper 10% critical values of | v |)* 


Upper 5% critical values of v= 


Ai 8} 


Ree 0-0 0-1 0-2 0-3 0-4 0-5 0-6 0-7 0-8 0.9 1.0 


f, 
6 104010 7-908 htl-860 1) 1-80") I:7e.N 1-740 A760 Leon 1-e58 1-01) 1-04 
8 1-94 9)01-60 11-85) | 1-80 | 1-76 0 673-4) 1-740) 1-760 1-79°1 1-92! 1 2-86 
10 | 1-94 | 1-90 | 1-85 | 1-80 | 1-76 | 1-73 | 1-73 | 1-74 | 1-76 | 1-78 | 1-81 
15 1-94 | 1-90 | 1-85 | 1-80 | 1-76 | 1-73 | 1-71 | 1-71 | 1-72 | 1-73 | 1-75 
20 1-G4 eld -O0a eal -S5. e180 | 1:76 1a 1-73 wip 17d lk 1-70n), 10a 1 alee 179 
© 1-94 | 1-90 | 1-85 | 1-80 | 1-76 | 1-72 | 1-69 | 1-67 | 1-66 | 1-65 | 1-64 
6 | 1:86 1-82.) 1-79 | 1-76 | 1-74 | 1-73 | 7-76 11-801 1-85°| 1-001 1-94 
8 1:86 | 1-82)! 1-70 | “1676 | 1-73 1 1-73 °4 1-78 | 1-78 | 1-79°7 1-62°7 2-86 
10 | 15869] 1-82) ) 179° } P1976 11-73 1-728 WP 72 7a) 1-760! 178!) 1-81 
15 1-86 | 1-82 | 1-79 | 1-76 | 1-73 | 1-71 | 1-71 | Le71 | 1-72 | 1-73 | 1-75 
20 1-86 | 1-82) | 1-79 | 1-76 | 1-73 | 1:714\) 1-700} 1-700) 1-700) 1-710) 4-92 
Ps 1-86 | |09-82)) 91-79) | 1-75 101-7249 I-700) 1-68 Wy 1-66Nl 1-655) 1-65) 1-64 
6 1810-780) 1-76: | Ie74 01°78) 1-730) 1-760ll 1804 1-850 1-90.) 1-04 
8 VSie i781 17601 1741 (70) 1-70 | 1730) i760) 1-700) soul 186 
10 1-$1)), 1-78.) 1-76 | 1-73 | 1-72 | W71| 1-72 | 173°) 1-76 4 1-78 || 1-81 
15 V80) | 1780) 1-76) 173° | 172 1 170 | 70 | 471 | 1-79] 1781 a-we 
20 1-61) 1:78" | 1-76 fiers fk F700) 160°) 1-69°") 1-70) 171) a2 
oo PSDP e7Sse 1-76 | 17S" | 171) 169-1 1-67") 1-66") 1-651 1-659) 1-62 
15 6 Lian 1-73) 1-72)) 17h | 7 |) 178 | 176 |} 180 | 2-85 | 1-00 | 1-04 
8 1-750 Seto? WV 1s7R 1-71 We Is7TL 0) L780 176") 2-700) 1-S20l) 1-86 
10 1750 Tonal Toul Tile TOR TOR 12er 1-7oh 176m) ive) el 
15 176A TON Al Tool 70 pe leTO ny L6dulT-70n) oe 1-7en a son ass 
20 1-75 los alt IC 170 e E9 4 69 4). 1-69.) 1-650) 1-G0a) Gail age 
00 1-75 | 1:73 | 112, | 1:70,|, 1-68 | 67.) VGbLl ese) 61650) eon a-Ge 
20 6 V2 aT Vo | 170° ' 1-70) 1-78" |" 1-765 1 1-80°) ese) 2-00) 1-04 
8 172 | 1-71 | 1-70 | 1:70 | 1-70.) Ievv | 1-73 | 1-76 | 1-79 | 1-82 | 1-86 
10 1-798 10-910 Pe1-70" |e1-GO ue 160 1-700) I-7t0! LTE F-7eml 17am 1-s8 
15 1790-71) 1701-60410 1-600) 1-6O0n 1-690) I-70 179i 12S) 1-78 
20 1-720) 1-712 el GO 160.0) 1-680 M 1-680) 1-68n) 160m) T7ON A-7d b-72 
00 1-72 | 1-71 | 1-70 | 1-68 | 1-67 | 1-66 | 1-66u) 65a 1Goe biped 
1-64 5 11-66. 1:67 |) 1 1:85 | 1-90 | 1-94 
1-64 5 | 166 | 66 | 1-6 1-79 | 1-82 | 1-86 
Bileites 1-660 1:76 | 1-78 | 1-81 
5 | 1-65 | 1-65 | 1. : 1-72 | 1-73 | 1-76 
Bova tees ei. ‘6 Lit Lan 1-72 
64 | 1-64 | 164 | 16 | Veal 364) 4-64 

| 
{ 


* y is normally distributed about 7 with variance A,o}+A,02, and s} and 6? are independent estimates 
of oj and a3, based on f, andf, degrees of freedom, respectively. A, and A, are known constants. _ 

In the problem of comparing the means of samples taken from two normal populations, put y= (%— 7)» 
f, = (n,—1),f,=(n,—-1), Ay=1/n, and A,= 1/ng, where n, and n, are the sample sizes. 


From Pearson and Hartley, (1966), pp. 142-145. 


oe ae 
: oe pas a a : 
_ : : - 

: ss. ~ ae 
| Cn, ™ 
oF ait. (ore EY Or work nl ‘Mos ene hiy errs oat 

Bitumen ylinaningae a Air 


Bi. ated =7 | wel we Twniten ae er a 
f } ’ y S - 
i °S i ' 


me it 1S. ehh! ibe por vt av ' ‘Tu oe i” i) 


————— ¥ ee ae as _ —— rs Sa = at 
. T. | ; 
Hl 810 Rail tf avoid ae | | pf taf } of | F.@ 
! } ; i] 
[eo eee eee | ee eee 
: - = 
ij . i . > | 
' i _ | 
\ gig i | i Mj j : 4 i i my -{ 
i isiG | ‘ i ’ ay! f ‘ i nl ‘ 
ue ey. | te | 
Le | . J H | ( | P ' & { 
i | i] : | 7 ! ‘ 
L') Ga i i i wt 
I } ; | 
¥ } 
a0.) i i a j 
f i 7 ' 
f { t + i 4 
t-] { j a vi 
i f eee | ‘ ' 5 poh 
+) ' Ti; a * 
re! i ) i 
{ ; 
Ly ee 
i] oes 7] 
I ei eat ae 
' la He H 
{ 5 mi | », oy rt vel aru 
f - | / L 
) 1 ney der d pee) ae) ep 
’ ‘ } 7 ? rad | 
' L ) ey BTM 
i ri wise) » % * 
oF | ce ae 71 mr) eed 
j nf j ' j j Pa | wn ime 
wl 
. : 
Toon is cP at) 6m | i 
ori Late | Of | yo) \ 
ee ea i iil iM Ort 1 et © 
} ar | i y, i) ; i f jaot 4 ns. t mis 
{ j OT i ; rer ) &e4 iovt 1 WEL | aps 
- F r . @ 
bist f ek | eg | 204) ai | OE VRE | be ee ae 
i ' 
ia} ra fa . iets ) { Ey | | Ay hel FE Gey | 
Ans | ea Gist Syt (al | & { sed Aha) 
1 <. ; ae ie eo 
’ | , “f cig uot ' Lise i vit f i-l } 
‘ } otf NT «] fil ry bled a amy © 
| wb. Se i . bay | 4 i bef ‘ay bay ( | i 
| | on 7 


—— ee a 


SAE 41) divi Prt) em @ ute bs | win Trai 7 fy 4 sanity ite 
. om! - yi » Peace tr » bron) LA. Aer to ‘{  yanel eat sin 
2s lg = : wT head ‘a Liezs Nelo Oved pln | @si.i wed Trae tee 


; ia wl 3 wen ie ira iy nt aH 


a . nad 


7. 
v 


8 


TableBl (continued) 


|| OR 0-1 0-2 0-3 0-4 


Asi | 
(Aya + A383) 


to bt ry 


f 

8 OS 1-25 no ON O14) e910) OS | 2-10 | 2-14 | 2-20 
10 31 |) 2-25 |) 2-80) S15 1-210) | 8-08) 208 1) Wo 
12 2-31 | 2:25 2-20 | 2-15 | 210 | 2:07 | 2-07 | 2-08 | 2-11 
15 O31 if 2-25) 1 2-20. 2-15. 1 2-10 | 2-07 4) 2-05 1 2-06 Il 2-08 
20 2:31 | 2:25 | 2:20. 12-15 | 62-10 152-06 | 2-04) 2-04 7 2-05 
co 2-31 | 2:25 | 2-20 | 2-14 | 2-09 | 2-05 | 2-01 | 1:99 | 1-97 
8 9-23 | 2-19 | 2-14 | 2:11 | 2:08 | 2-08 | 2-10 | 2-15 | 2-20 
10 2-23 | 2:18 | 2-24 | 2-11 | 2-08 | 2-06 | 2-08 | 2-11 | 2-14 
2 993° [2418 | 214 |°2-10 12-07 | 2-06 | 2-06 | 2-08 | 2-11 
15 9.23 | 2-18 | 2-14 | 2:10 | 2-07 | 2:05 | 2-05 | 2-06 | 2-08 
20 9S) O16) (e940 189-01 2:07 1) (2:05 12-04 I Odo 2-05 
© O23 | O18 12:14, | 2-10 | “2-06 | 2:03 | 2-00 | 1-98 1-97 
8 FIG 1 2s ip 2-2 (22-0811 207 | 07 | 2-10 | 2:15 | 2-20 
10 2-18 | 2-14 | 2-11 | 2-08 | 2:06 | 2-06 | 2-07 | 2:10 | 2-14 
12 9-18 | 2-14 | 2:11 | 2-08 | 2:06 | 2-05 | 2-06 | 2-08 | 2-11 
15 9-18 | 2:14 12-11 |°2-08 | 2-06 | 2-04 | 2-04 | 2-06 | 2-08 
20 2-18 | 2-14 | 2-11 | 2-08 | 2-05 | 2-04 | 2-03 | 2-03 | 2-05 
ee) 2-18 | 2-14 | 2-11 | 2-07 | 2-04 | 2-02 | 1:99 | 1-98 | 1-97 
8 2-13 | 2-10 | 2-08 | 2-06 | 2:05 | 2-07 -| 2:10 | 215 | 2-20 
10 2:13 | 2:10 | 2-08 | 2-06 | 2-05 | 2-05 | 2-07 | 2-10 | 2:14 
12 2-13 | 2-10 | 2-08 | 2-06 | 2-04 | 2-04 | 2-06 | 2-08 | 2-11 
15 2-13 | 2:10 | 2:08 | 2:05 | 2-04 | 2:03 | 2:04 | 2-05 | 2-08 
20 2-13 | 2-10 | 2-08 | 2-05 | 2:04 | 2-03 | 2-03 | 2-03 | 2-05 
re 2-13 | 2:10 | 2-07 | 2:05 | 2:02 | 2-00 -99 | 1:97 | 1-97 
8 9:09 | 2-07 | 2-05 | 2-04 | 2-04 | 2-06 | 2-10 | 9-15 | 2-20 
10 209 | 2-06 | 2-05 | 2-04 | 2:04 | 2-05 | 2:07 | 2:10 | 2-14 
12 2-09 | 2-06 | 2-05 | 2-03 | 2-03 | 2-04 | 2-05 | 2-08 | 2-11 
15 2-09 | 2:06 | 2:05 | 2-03 | 2-03 | 2-03 | 2-04 | 2-05 | 2-08 
2 2:09 | 2-06 | 2:05 | 2-038 | 2-02 | 2:02 | 2:02 | 2:03 | 2-05 
oO 2:09 | 2-06 504.01 52502 0153-050) 2 1:98 | 1:97 | 1:96 
8 1:96 | 1-96 | 1-97 | 1-99 | 2-01 | 2-05 | 2-09 | 2-14 | 2-20 
10 1-96 | 1-95 | 1-97 | 1-98 | 2-00 | 2-03 | 2-06 | 2-10 | 2-14 
12 1-96 | 196 | 1-97 | 1-98 | 1:99 | 2-02 | 2-04 | 2-07 | 2-11 
15 1-96 | 1:96 | 1-97 | 1:97 | 1-99 | 2-00 | 2-02 | 2-05 | 2-07 
20 1-96 | 1-96 | 1-96 | 1-97 | 1-98 | 1-99 | 2-01 | 2-02 | 2-04 
oo 1:96 | 1-96 | 1-96 | 1-96 | 1:96 | 1:96 | 1:96 | 1-96 | 1-96 
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TableBl (continued). Test for comparisons involving two variances which must be 
separately estimated 


(y-7) 
(Ay 8} + A283) 
(i.e. upper 2% critical values of | v |)* 


Upper 1% critical values of v= 


2 
AG os 


Asi t+ A285 


ae 

10 2-76 | 2-70 | 2-63 | 2-56 | 2-51 { 2-50 | 2-51 | 2-56 2:70 | 2-76 
2-76 | 2-70 | 2-63 | 2-56 | 2-51 | 2-49 | 2-49 | 2-52 2-62 | 2-68 

2:76 | 2:70 | 2-63 | 2-56 | 2-51 | 2-48 | 2-47 | 2-48 2-56 | 2-60 

2-76 | 2-70 | 2-63 | 2-56 | 2-51 | 2-47 | 2-45 | 2-45 2:49 | 2-53 

2:76 | 2-70 | 2-63 | 2-56 | 2-50 | 2-46 | 2-43 | 2-42 2-44 | 2-46 

2-76 | 2-70 | 2-63 | 2-56 | 2-50 | 2-44 | 2-40 | 2-36 2-33 | 2-33 

12 2-68 | 9-62 | 2-57 | 2:52 | 2-49 | 2-49 | 2-51 | 2-56 2-70 | 2-76 
2-68 | 2-62 | 2-57 | 2-52 | 248 | 2-47 | 2-48 | 2-52 2:62 | 2-68 

| 9-68 | 2-62 | 2:57 | 2-52 | 2-48 | 2-46 | 2-46 | 2-48 | 2-52 | 2-56 | 2-60 

2-68 | 2-62 | 2-57 | 2-52 | 2-48 | 2-45 | 2-44 | 2-45 | 2:47 | 2-49 | 2-53 

2-68 | 2-62 | 2-57 | 2-52 | 2-47 | 2-44 | 2-42 | 2-41 | 2-42 | 2-44 | 2-46 

2-68 | 2-62 | 2:67 | 2-51 | 2-46 | 2-42 | 2-38 | 2-36 | 2-34 | 2-33 | 2-33 

15 10 2-60 | 2:56 | 2-52 | 2-48 | 2:47 | 2-48 | 2-51 | 2-56 | 2-63 | 2-70 | 2-76 
12 2-60 | 2-56 | 2-52 | 2-48 | 2-46 | 2-46 | 2-48 | 2-52 | 2-57 | 2-62 | 2-68 

15 2-60 | 2:56 | 2-51 | 2-48 | 2-45 | 2-45 | 2-45 | 2-48 | 2-61 | 2-56 | 2-60 

20 2-60 | 2:56 | 2-51 | 2-48 | 2-45 | 2-43 | 2-43 | 2-44 | 2-46 | 2-49 | 2-53 

30 2:60 | 2-56 | 2-51 | 2-47 | 2:44 | 2-42 | 9-41 | 2-41 | 2-42 | 9-44 | 2-46 

© 9-60 ; 2°56 | 2-51 | 2:47 | 2-43 | 2-40 | 2-37 | 2-35 | 2-34 | 2-33 | 2-33 

20 10 9-53 | 2-49 | 2-47 | 2-45 | 2-45 | 2-47 | 2-51 | 2-56 | 2-63 | 2-70 | 2-76 
12 2:53 | 2-49 | 2-47 | 2-45 | 2-44 | 2-45 | 2-48 | 2-52 | 2-57 | 2-62 | 2-68 

15 9:53 | 2-49 | 2-46 | 2-44 | 2-48 | 2-45 | 2-45 | &-48 | 2-51 | 2-56 | 2-60 

20 | 2-53 | 2-49 | 2-46 | 2-44 | 2-42 | 2-42 | 2-42 | 2:44 | 2-46 | 2-49 | 2-53 

30 | 9-53 | 2:49 | 2-46 | 2-44 | 2-42 | 2-40 | 2-40 | 2-40 | 2-42 | 2-43 | 2-46 

© 2-53 | 2:49 | 2-46 | 2-43 | 2-40 | 2-38 | 2-36 | 2-34 | 2:33 | 2-33 | 9-33. 

30 10 2:46 | 2-44 | 2-42 | 2-42 | 2-43 | 2-46 | 2-50 | 2-56 | 2-63 | 2-70 | 2-76 
12 2-46 | 2-44 | 2-42 | 2-41 | 2:42 | 2-44 | 2-47 | 2-59 | 2-57 | 2-62 | 2-68 

15 9:46 | 9:44 | 2:42 | 9-41 | 2-4] | 2-42 | 2-44 | 2-47 | 9-51 | 2-56 | 2-60 

20 2-46 | 2:43 | 2-42 | 2-40 | 2-40 | 2-40 | 2-42 | 2:44 | 2-46 | 2-49 | 2-53 

30 2:46 | 2-43 | 2:42 | 2-40 | 2-39 | 2:39 | 2-39 | 2-40 | 2-42 | 2-43 | 2-46 

ree) 2-46 | 2:43 | 2:41 | 2:39 | 2:37 | 2-36 | 2-35 | 2-34 | 2:33 | 2:33 | 2-33 

co 10 2-33 | 2-33 | 2-34 | 2-36 | 2-40 | 2-44 | 2-50 | 2-56 | 2-63 | 2-70 | 2-76 
12 2-33 | 2-33 | 2:34 | 2-36 | 2-38 | 2-42 | 2-46 | 2-51 | 2-57 | 2-62 | 2-68 

15 9-33 | 2-38 | 2-34 | 2-35 | 2-37 | 2-40 | 9-43 | 9-47 | 2-51 | 2-56 | 2-60 

20 2-33 | 2-33 | 2:33 | 2-34 | 2-36 | 2-38 | 2-40 | 2-43 | 2-46 | 2-49 | 2-53 

30 O38 |) Opa) 283. 1) 294") 2-35. |) 8-36 | 2:37 | 2:39 | 2-41 | 2-43 | 2-46 
ee) 2-33 | 2-33 | 2-33 | 2-33 | 2-33 | 2-33 | 2-33 | 2.33 | 2-33 | 2-33 | 2-33 | 


* y is normally distributed about 7 with variance A,o}+A,03, and sj and s3 are independent estimates 
of oj and o3, based on f, and f, degrees of freedom, respectively. A, and A, are known constants. 

In the problem of comparing the means of samples taken from two normal populations, put y=(Z%,—2;,), 
f,=(n,—1), £=(n,—-1), Ay=1/ny and A,=1/n,, where n, and n, are the sample sizes. 
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Table B1 (continued) 


Upper 0-5 % critical values of v= - he IE 
/o f A AgS: + A, 83) 


(i.e. upper 1% critical values of | v |) 


2 
ee 0-1 1-0 
(Ay s} + A283) 
CRU ee ea et ae oe eee Se 
| | ! 
f ae! | 
10 3-17 | 3-08 ; 3-00 | 2-90 | 2-82 | 2:79 ; 2-82 | 2-90 | 3-00 | 3-08 | 3-17 
12 3:17 | 3-08 | 3-00 | 2-91 | 2-82 | 2-78 | 2-79 | 2-84 | 2-91 | 2-98 | 3-05 
15 3-17 | 3-08. | 3-00 | 2-91 | 2:82 | 2-77 | 2-76 | 2-78 | 2-83 | 2-89 | 2-95 | 
20 3:17 | 3:08 | 3-00 | 2-91 | 2-82 | 2-76 | 2-73 | 2-74 | 2-76 | 2-80 | 2-85 
30 3-17 | 3-08 | 3-00 | 2-91 ; 2-82 | 2-75 | 2-71 | 2-69 | 2-70 | 2-72 | 2-75 
oO 3-17 | 3:08 | 2-99 | 2-91 | 2-82 | 2-74 | 2-67 | 2-63 | 2:60 | 2-58 2-58 | 
10 3-05 | 2-98 | 2-91 | 2-84 | 2-79 | 2-78 | 2-82 | 2-91 | 3-00 | 3-08 | 3-17 
12 3:05 | 2:98 | 2:91 | 2-84 | 2:78 | 2:76 | 2-78 | 2-84 | 2-91 | 2:98 | 3-05 
15 3-05 | 2-98 | 2-91 | 2-84 | 2-78 | 2-75 | 2-75 | 2:78 | 2-83 | 2-89 | 2-95 | 
20 3-05 | 2:98 | 291 | 2-84 | 2-78 | 2-74 | 2-72 | 2-73 | 2-76 | 2-80 | 2-85 
30 3:05 | 2:98 | 2:91 | 2-84 | 2:77 | 2-73 | 2-70 | 2-69 | 2-70 | 2-72 | 2-75 
oO 3-05 | 2-98 | 2:91 | 2-84 | 2-77 | 2-71 | 2-65 | 2:62 | 2:59 | 2:58 | 2-58 
10 2-95 | 2-89 | 2-83 | 2-78 | 2-76 | 2-77 | 2-82 | 2-91 | 3-00 | 3-08 | 3-17 
12 2-95 | 2-89 | 2-83 | 2-78 | 2-75 | 2-75 | 2-78 | 2-84 | 2-91 | 2-98 | 3-05 
15 2-95 | 2-89 | 2-83 | 2-78 | 2-74 | 2-73 | 2-74 | 2-78 | 2-83 | 2-89 | 2-95 
20 2-95 | 2-89 | 2-83 | 2-78 274) O27 | 2734) 2-76 | 2-80 | 2-85 
30 2-95 (-2-80) 1) 2°62) 1 62°78) | 2-73.) 2°70 1, 2-68 12-68 152-70 | 12-720 9-75 
0 2-95 | 2-89 | 2-83 | 2-77 | 2-72 | 2:67 | 2-64 | 2:61 | 2:59 | 2-58 | 2-58 
10 2-85 | 2:80 | 2-76 | 2-74 | 2-73 | 2-76 | 2-82 | 2-91 | 3-00 | 3-08 | 3-17 
12 2-85 | 2-80 | 2-76 | 2:73 | 2-72 | 2-74 | 2-78 | 2-84 | 2-91 | 2-98 | 3-05 
15 9-85 | 2:80 | 2:76 | 2-73 | 2-71 | 2-71 | 2-74 | 2-7 2:83 | 2:89 | 2-95 
26 2°85.) 2-80 | 2:76 192-73 | 2-70.) 2-70 | 2-70) 2-73.) 2-76 1 2-80) | 2-85 
30 2:85 || 2-80 | 2-761 2-72 | 2:60 | 2-68 | 2-67 | 2-68 | 2:70 | 2-72 | 2:75 
roe) 2-85 | 2-80 | 2-76 | 2-72 | 2-68 | 2-65 | 2-62 | 2-60 | 2-59 | 2-58 | 2-58 
10 2-75 | 2-72 | 2-70 | 2-69 | 2-71 | 2:75 | 2-82 | 2-91 | 3-00 | 3-08 | 3-17 
12 2:75 | 2-72 | 2-70 | 2-69 | 2-70 | 2-73 | 2:77 | 2-84 | 2-91 | 2-98 | 3-05 
15 2:75 | 2-72 | 2-70 | 2-68 | 2-68 | 2-70 | 2:73 | 2:78 | 2-83 | 2-89 | 2-95 
20 2-15 | 2-72 | 2-70 | 2°68 | 2-67 | 2-68 | 2-69 | 2-72 | 2-76 | 2-80 | 2-85 
30 2-75 | 2-72 | 2-69 | 2-67 | 2-66 | 2-66 | 2-66 | 2-67 | 2-69 | 2-72 | 2-75 
oo 2:75 | 2:72 | 2-69 | 2-66 | 2-64 ! 9-62 | 2:60 | 2:59 | 2-58 | 2-58 | 2:58 | 
10 2-58 | 2:58 | 2-60 | 2-63 | 2-67 | 2-74 | 2-82 |. 2-91 | 2-99 | 3-08 | 3-17 
12 2:58 | 2-58 | 2-59 | 2-62 | 2-65 | 2-71 | 2-77 | 2-84 | 2-91 | 2-98 | 3-05 
15 2-58 | 2-58 | 2-59 | 2-61 | 2-64 | 2-67 | 2:72 , 2:77 | 2-83 | 2-89 | 2-95 
20 2-58 | 2-58 | 2-59 | 2-60 | 2-62 | 2-65 | 2-68 | 2-72 | 2:76 | 2-80 | 2-85 | 
30 2-58 | 2-58 | 2-538 | 2-59 | 2:60 | 2-62 | 2:64 | 2:66 | 2-69 | 2-72 | 2-75 
©O 2:58 | 2-58 | 2-58 | 2-58 | 2-58 | 2-58 | 2-58 | 2-58 | 2-58 | 2-58 | 2-58 
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